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Abstract 



New measurements are presented of the photon structure function F'^ix^Q'^^ at four 
values of between 9 and 59 GeV^ based on data collected with the OPAL detec- 
tor at centre-of-mass energies of 161—172 GeV, with a total integrated luminosity of 
18.1 pb~^. The evolution of with in bins of x is determined in the range 
from 1.86 to 135 GeV^ using data taken at centre-of-mass energies of 91 GeV and 
161 — 172 GeV. F^ is observed to increase with with a slope of a~^dF^/dlnQ^ = 
O.IOT 

ao3 measured in the range 0.1 < x < 0.6. 
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1 Introduction 



The measurement of the photon structure function F^, and in particular of its evolution 
with the momentum transfer squared, Q^, is a classic test of perturbative QCD 
is expected to increase only logarithmically with Q^, therefore the large range of 
values accessible at the e"'"e~ collider LEP, which presently extends from about 1 GeV^ 
to about 400 GeV^ and which will increase in future, makes it an ideal place to study 
this evolution. 

The results presented in this paper are based on new measurements of F2 using 
data in the range from 6 to 100 GeV^ recorded by the OPAL detector in 1996 at 
y/s^ = 161 — 172 GeV and on results of the analysis of data collected at e+e~ 
centre-of-mass energies ^/s^ around the mass of the Z° (denoted by ^/s^ = 91 GeV). 
This measurement is an extension of the analysis of detailed in using basically 
the same methods to analyse the singly-tagged two-photon events at higher ^/s^. 

For singly-tagged event^ the process e"'"e~ — > e"'"e~ -|- hadrons can be regarded 
as deep inelastic scattering of an e^ on a quasi-real photon, where the flux of quasi- 
real photons can be calculated using the equivalent photon approximation The 
cross-section for deep inelastic electron-photon scattering is expressed as 

= 1^ [(1 + (1 - yf) Flix, - y^Flix, Q')] (1) 

where = —q^ is the negative value of the four-momentum squared of the virtual 
photon. The usual dimensionless variables of deep inelastic scattering, x and y, are 
defined as x = Q^/2(p ■ q) and y = {p ■ q)/{p ■ k), and a is the fine structure constant. 
The symbols p, q and k denote the four-vectors of the quasi-real photon, the virtual 
photon and the incoming electron respectively. The structure function F^ is related to 
the sum over the parton densities of the photon (see e.g. f^). In the kinematic region 
of low values of y studied {y"^ -C 1) the contribution of the term proportional to the 
longitudinal structure function Fl(x, Q^) is negligible. 

The paper is organised as follows. After the description of the OPAL detector in 
section Q the data selection and background estimates are detailed in section |^. The 
quality of the description of the observed hadronic final state by the Monte Carlo 
models is discussed in section ^. The measurement of the evolution of F2 is outlined 
in section ||. The results are presented in section |, followed by the conclusions given 
in section |^ 



^The term singly-tagged denotes the situation where the electron, which radiates the virtual photon, 
undergoes deep inelastic scattering and is seen (tagged) in the detector, whereas the other electron, 
which radiates the quasi-real photon, is unseen as it leaves the detector close to the beam direction. 
In this paper positrons are also referred to as electrons, and the electron and positron masses are 
neglected. 
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2 The OPAL detector 



The OPAL detector is described in detail elsewhere 0; only the subdetectors which 
are most relevant for this analysis, namely the electromagnetic calorimeters and the 
tracking devices, are detailed below. The OPAL detector has a uniform magnetic 
field of 0.435 T along the beam direction throughout the central tracking region, with 
electromagnetic and hadronic calorimetry and muon chambers outside the coil. 

In the OPAL right-handed coordinate system the x-axis points towards the centre 
of the LEP ring, the y-axis points upwards and the 2;-axis points in the direction of the 
electron beam. The polar angle 6 and the azimuthal angle are defined with respect 
to the 2-axis and respectively. 

The small-angle silicon tungsten luminometer (SW) covers the region in 6 from 
25 to 59 mrad. From 1996 onwards, the lower boundary of the acceptance has been 
increased to 32 mrad following the installation of a low angle shield to protect the 
central detector against possible synchrotron radiation. The SW consists of two finely 
segmented calorimeters placed around the beam pipe, one each side of the detector, 
at a distance of 2.4 m from the interaction point. Each calorimeter is composed of a 
stack of 19 silicon wafers interleaved with 18 tungsten plates placed perpendicular to 
the beam axis. The sensitive area of each silicon wafer extends in radius from 62 to 
142 mm from the beam axis and is segmented into 32 pads radially and 32 in azimuth 
around the beam. The radial position of electron showers in the calorimeter can be 
determined with a typical resolution of 0.06 mrad in 6. The energy resolution for a 
beam energy electron of 86 GeV energy is approximately 2.5 GeV. 

The clean acceptance of the forward detectors (FD) covers the 6 region from 60 
to 140 mrad at each end of the OPAL detector. The FD consists of cylindrical lead- 
scintillator calorimeters with a depth of 24 radiation lengths {Xq) divided azimuthally 
into 16 segments. The energy resolution for electromagnetic showers is 18%/\/E, where 
E is in GeV. An array of three planes of proportional tubes buried in the calorimeter 
at a depth of 4 Xq provides a precise shower position measurement, with a typical 
resolution of 3-4 mm, corresponding to 2.5 mrad in 6, and no more than 3.5 mrad in 
0. 

Charged particles are detected by a silicon microvertex detector, a drift chamber 
vertex detector, and a jet chamber. Outside the jet chamber, but still in the magnetic 
field, lies a layer of drift chambers whose purpose is to improve the track reconstruction 
in the z-coordinate. The resolution of the transverse momentum for charged particles 
is cTpJpt = a/(0.02)2 _|_ (o.0015pt)^ for | cos6'| < 0.7, where pt is in GeV, and degrades 
for higher values of | cos 6 \ . 

Both ends of the OPAL detector are equipped with electromagnetic endcap calorime- 
ters covering the range from 200 to 630 mrad in polar angle. They are homogeneous 
devices composed of arrays of lead-glass blocks of 9.2 x 9.2 cm^ cross-section and typ- 
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ically 22 Xq in depth, giving good shower containment. In the central region, outside 
the solenoid, is the electromagnetic barrel calorimeter of similar construction. 

The deep inelastic scattering events are triggered with high efficiency by the large 
energy deposits of the scattered electron in the calorimeters and by charged particle 
tracks seen in the tracking devices. 



3 Kinematics and data selection 

The measurement of i^2^(x, Q"^) requires the determination of x and which are ob- 
tained from the four-vectors of the tagged electron and the hadronic final state as 
follows: 

Q2 = 2EbStag(l-COSetag) (2) 



Here E'tag and 6'tag are the energy and polar angle of the tagged electron, £'b is the 
beam energy, and W the invariant mass of the hadronic final state. = —p^ is the 
negative value of the virtuality of the quasi-real photon. For this singly-tagged sample, 
an antitag condition on the second electron is applied (see list of cuts below). This 
ensures that P^ is much smaller than and it is therefore neglected in the evaluation 
of X using Eq. ^ 

The four-momentum of the hadronic system is calculated by summing over all 
charged particle tracks (assuming the pion mass) and calorimeter clusters, where qual- 
ity criteria are applied to both the tracks and the clusters to ensure that they are well 
reconstructed. To avoid double counting of particle momenta, a matching algorithm 
for the association of a charged particle track with a cluster is applied 0]. 

The analysis uses the data at ^/s^ = 161 — 172 GeV with an integrated e"'"e~ 
luminosity of 18.1±0.1 pb^^, as determined from small-angle Bhabha scattering events. 
The tagged electron is detected either in the SW detectors (Q^ ~ 6 — 20 GeV^) or in 
the FD detectors (Q^ f» 20 — 100 GeV^). These two samples are subject to slightly 
different selection criteria and are referred to as the SW and FD samples. Candidate 
events for the process 7*7 — > hadrons are required to satisfy criteria for the tagged 
electron as well as for the hadronic final state, in addition to several technical cuts, 
to ensure good detector status and track quality. The event selection listed below is 
designed to have a high efficiency for signal events and to reject background events, 
which mainly stem from the process 7*7 — >■ t^t~ . The requirements for the SW (FD) 
samples are: 

1. A tagged electron candidate is required to produce a cluster in a SW (FD) de- 
tector with an energy of E'tag > 0.775 i?b (-f'tag > 0.60 i?b) and a polar angle 
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of 33 < 6'tag < 55 (60 < 6'tag < 120) mrad with respect to either of the beam 
directions. 

2. The energy E'at of the most energetic cluster in the hemisphere opposite to the 
one which contains the tagged electron is restricted to -Eat < 0.08 -Eb (antitag 
requirement). 

3. There must be at least three tracks originating from the hadronic final state, 
iVch > 3. 

4. The visible invariant mass W^is of the hadronic system, calculated as the mass of 
the four-momentum vector of the hadronic system as defined above, is required 
to be in the range 2.5 < Wvis < 40 GeV. 

The harder requirement on the tag energy applied to events of the SW sample reflects 
the much higher background from off-momentum electrons closer to the beam direction. 
An off-momentum electron is an electron that was lost from the beam and scattered 
into the detector. In addition to the off-momentum electron, due to an interaction that 
occurs close to the nominal interaction point, charged particles are produced which 
fulfill the requirements for the hadronic final state. Since the energy spectrum of the 
off-momentum electrons peaks around 0.5i?b this background is essentially eliminated 
by requiring Et^g > 0.775 -Eb- 

With these cuts 879 and 414 events with average squared momentum transfers 
(Q^) of approximately 11 GeV^ and 41 GeV^, are selected in the SW and FD samples 
respectively. The accessible x range for the two samples is 0.004 < x < 0.76 and 
0.012 < X < 0.94 respectively. The trigger efficiency is evaluated from the data to be 
above 98% for the SW sample and essentially 100% for the FD sample. 

The background to the 7*7 — hadrons signal comes from events which contain a 
true or fake tagged electron and an apparent low-mass hadronic final state (compared 
to y^i^). The dominant source of this background is 7*7 — > leptons and particularly 
7*7 T^T^. These processes were simulated with the Vermaseren program [Q. The 
background contribution from 7*7 — ^ t~^t~ events is 36.0 ± 1.8 (29.0 ± 2.2) events 
and the 7*7 — *• e"'"e~ process contributes 16.8 it 1.2 (9.1 ± 0.9) events to the SW (FD) 
samples respectively. Additional background sources like hadrons events, four- 

fermion events with e"'"e~ qq and e"'"e~ i~^£~ {£ = e, fi, r) final states, and Z° t^t^ 
events, were studied. The total contribution of these processes to the background is of 
the order of 1 and 3 events in the two samples, respectively, and consequently they are 
neglected. 

Events with a scattered electron observed in the electromagnetic endcap calorimeter 
allow higher values of to be reached than in the SW and FD samples. With the 
current level of integrated luminosity collected at ^/s^ =161 — 172 GeV too few such 
events have been observed to allow a detailed analysis to be performed. 
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4 Description of the hadronic final state 



The measurement of requires the determination of W from the hadronic final state. 
Because of the finite detector resolution and the incomplete angular coverage, especially 
in the forward region, the correlation between visible hadronic mass, Wyis, and W 
depends critically on the modelling of the hadronic final state, which is correlated with 
the observed hadronic energy flow. Therefore a detailed comparison of the observed 
hadronic final state and the predictions from the various Monte Carlo models is needed. 
The results of this study, which follows closely the one described in ^ , are summarised 
in this section. 



The set of Monte Carlo generators used to simulate the signal events differs only 
slightly from that used in the previous analysis 0. HERWIG |^ version 5.9 is used 
instead of 5.8d with the charm quark mass altered from 1.8 GeV to 1.55 GeV. All 
events are generated without the simulation of the soft underlying event 0. The 
F2GEN [|10| generator is used to simulate events based on F2 for four flavours (u,d,s,c) 
with rric = 1.5 GeV, rather than with three flavours and with the charm contribution 
added using the Vermaseren program, as in . The PYTHIA [|1^ Monte Carlo remains 
the same as in . The integrated luminosity of each Monte Carlo sample corresponds 
to approximately ten times the integrated luminosity of the data. The Monte Carlo 
generation was based either on the leading order (LO) GRV [0 or the LO SaSlD [l^ 
parton density parametrisations. All Monte Carlo events are passed through the OPAL 
detector simulation program and the same reconstruction and analysis chain as 
the real data events. 



For the comparison of the observed hadronic final states to the Monte Carlo pre- 
dictions a particle is defined as a track or a cluster resulting from the matching al- 
gorithm 1^. In PI it was shown that none of the models used accurately described 
the angular distribution and the transverse energy of the hadronic final state. These 
discrepancies persist at higher centre-of-mass energies. The accessible region for the 
SW sample is similar to that of the low-Q^ sample of the analysis at ^/s^ = 91 GeV p|, 
where the scattered electron was observed in the FD. 

Figure |l| shows the hadronic energy flow per event, 1/A^ dE/drj, as a function 
of pseudorapidity, rj = — ln(tan(6''/2)), for the SW sample and the FD sample for 
two bins in Xvis- Only statistical errors are shown. The tagged electron is always at 
negative pseudorapidity and is not shown. In the region 0.1 < Xvis < 0.6 the agreement 
between data and the various models is satisfactory, with the exception of the region 
—4 < 7] < —2 for the SW sample where all Monte Carlos are low compared to the 
data. This effect comes from leakage into the FD detector and is almost entirely 
eliminated if the requirement for the maximum 6'tag for the SW sample is restricted 
to smaller angles. This is taken into account in the evaluation of the systematic error 
on F2 . Significant differences are observed for x^is < 0.1. In the remnant direction at 
positive T], both HERWIG and PYTHIA tend to overestimate the energy deposited in 
the forward region of the detector. In the central region, however, they underestimate 
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the transverse energy. In contrast, the pointhke events generated with the F2GEN 
model he higher than the data in the central region and lower in the forward region. 
The number of events with Xvis > 0.6 is insufficient to allow firm conclusions to be 
drawn. 

The minimum tag energy and the antitag requirement for the FD sample of this 
analysis are slightly different from the ones used for the low-Q^ sample in This 
results in a different shape of the energy fiow distributions and a slightly improved 
quality of the data description by the Monte Carlo models. 

In summary, in the intermediate range of x^is a good overall agreement between 
the data and the Monte Carlo models is found for both samples. At x^is < 0.1 the 
fiow of hadronic energy into the central region of the detector is underestimated by 
the HERWIG and PYTHIA Monte Carlo models, and overestimated by the F2GEN 
model. Improved Monte Carlo models would be therefore required in order to reduce 
the model-dependent systematic errors on the measurement of in this region. 



5 Determination of F2 

To obtain the photon structure function in bins of from the measured Xvis 
distribution the method of regularised unfolding is employed |jl5|. In this method. 



a reference Monte Carlo is used to obtain the detector efficiency and resolution and 
the regularisation reduces the effect of inherent oscillating fiuctuations of the unfolded 
structure function. See for more details of this procedure. To unfold F2 from the 
data HERWIG, with the structure function based on the GRV parametrisation, 
is chosen as the reference Monte Carlo. This choice is motivated by the fact that 
HERWIG gives a slightly better description of the hadronic final state than PYTHIA, 
(see and Fig. |l]). The reference Monte Carlo is used to determine the central values 
and most of the systematic uncertainties. The PYTHIA and the F2GEN models are 
used to estimate the model dependence of the result. 

The data at ^/s^ = 161 — 172 GeV are subdivided into four ranges, two ranges 
for the SW sample and two for the FD sample. The ranges in are 6-11, 11-20, 20- 
40 and 40-100 GeV^ with average squared momentum transfers of (Q^) = 9, 14.5, 30 
and 59 GeV^ respectively. The Xvis distributions of the SW and FD samples are shown 
in Figs. I^a and 0b. The distributions of the reference Monte Carlo with background 
added are shown by the dashed histograms. In addition, the background events are 
shown separately at the bottom of the figure. The number of Monte Carlo events 
is the absolute prediction for the data luminosity. The solid histogram represents 
the reweighted distribution of the signal Monte Carlo with background added after 
the unfolding has been performed 0. From the figures one observes that the mean 
Xvis increases with increasing Q^, and that the a;vis distribution of the data is well 
represented by the sum of the reweighted signal Monte Carlo and background Monte 
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Carlo samples after unfolding. 

The weight factors for each individual event, obtained from the unfolding, can be 
used to construct reweighted Monte Carlo distributions of different variables. The 
reweighted distributions are the Monte Carlo predictions based on the structure func- 
tion unfolded from the data. Any reweighting based on the generated x distribution 
will change the shape of other measurable variables besides Xvis- As all variables are 
differently correlated with x^is the comparison of their distributions with the data gives 
an important check of the transformation, as described by the Monte Carlo simulation, 
between the partonic distributions and the measurable distributions. As an example 
Figs. ^ and ^ show the distributions of the SW and FD samples, applying the 
cuts listed in Sect. ^. There is good agreement in the shape of the distributions. The 
distributions from both Monte Carlo samples differ slightly in normalisation from the 
distributions observed in the data, and the agreement improves after the reweighting 
based on the unfolding of the Xvis distribution. A similar behaviour is found for other 
observables, both for quantities measured from the electron, like i^tag and ^tag, and for 
distributions obtained from the hadronic final state, like N^h and VTvis- 



6 Results 



The unfolded F2 measurements for the data taken at y/s^ = 161 — 172 GeV are shown 
in Figs. 0a— and are listed in Table |I|. The value of F2 /a is shown at the centre 
of the X bin. The bin sizes are indicated by the vertical lines at the top of the figure. 
The error bars show both the statistical error alone and the full error, given by the 
quadratic sum of statistical and systematic errors. The central values and statistical 
errors of the F2 measurements are estimated by using the reference unfolding. The 
F2 values presented here are not corrected for the fact that is not strictly equal to 
zero 0. 

The estimation of the systematic error includes three parts : the variation of the 
compositions of signal and background events in the sample, the use of different 
structure functions assumed in the Monte Carlo samples, and the different modelling 
of the formation of the hadronic final state. To allow for a different composition of 
signal and background events, the event selection cuts are varied. The choice of the 
cut variations reflects the different population of signal events in the four ranges, in 
terms of the scattering angle of the electron and W^is, as well as the different behaviour 
of the background events. In each case the unfolding is carried out using the reference 
Monte Carlo model where only one cut is varied from the standard set. To study the 
uncertainty due to the structure functions assumed in the Monte Carlo samples, for 
the SW sample, the unfolding is done using the HERWIG generator, the standard set 
of cuts and the SaSlD parton density parametrisations. The effect of the different 
modelling of the formation of the hadronic final state is studied by repeating the 
unfolding using the PYTHIA Monte Carlo model with the SaSlD parametrisation, 
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and the F2GEN Monte Carlo model with the GRV parametrisation. The systematic 
error assigned to the result, shown in Table |l], is taken as the maximum deviation 
of any unfolding result from the central values, which in some cases leads to rather 
asymmetric errors, particularly when the HERWIG Monte Carlo model is replaced by 
PYTHIA or F2GEN. The systematic error is dominated by the model uncertainties. 

The measured F2 as function of x is almost flat within the region studied and 
the absolute normalisation of F2 is well described by various predictions. Shown in 
Fig. 1^ are the F2 predictions for the leading order GRV and the SaSlD parton den- 
sity parametrisations evaluated at the corresponding values of (Q^). In each case the 
expected contribution from massive charm quarks is added. The fact that the charm 
threshold depends on x and varies with causes the change in shape of the pre- 
dictions, see Eq. || for 14^ = 2mc. Both parametrisations successfully describe the 
measured function of x in all Q ranges. 

Figure |^ also shows an augmented asymptotic prediction for F^. The contribution 
to F2 from the three light flavours is approximated by Witten's leading order asymp- 
totic form [|l|, using the parametrisation given in [0. This has been augmented by 
adding a charm contribution evaluated from the Bethe-Heitler formula |T^, and an 



estimate of the hadronic part of the photon structure function, which essentially cor- 
responds to the hadronic part of the GRV (LO) parametrisation of F^. The hadronic 



contribution is derived from the structure function parametrisation of the pion 



using VDM arguments as detailed in and evolved to the corresponding values of 
(Q^). The components are evaluated in leading order using A^^ = 0.232 GeV. It is 
known that the asymptotic solution, which is obtained by neglecting terms in 

the evolution equations that vanish for Q'^ ^ 00, has deficits in the region of low x, be- 
cause the remaining terms create divergences in the solution which do not occur in the 
solution of the full evolution equations ||T9[. However, the asymptotic solution has the 



appealing feature that it is calculable in QCD, even at higher order and for medium 
X and with increasing it should be more reliable. In addition at high x and the 
hadronic contribution is expected to be small. In the region of medium x values studied 
here this asymptotic prediction in general lies higher than the GRV and SaS predictions 
but it is still in agreement with the data. The importance of the hadronic contribution 
to Fj, which is shown separately at the bottom of the figure, decreases with increasing 
X and Q^, and it accounts for only 15% of F2 at = 59 GeV^ and x = 0.5. The 
asymptotic solution increases with decreasing A|^^. For = 59 GeV^ and x = 0.5 the 
change in F^ is +24% and -16% if Af^ is changed from Af^ = 0.232 GeV to 0.1 GeV 
and 0.4 GeV respectively. 

As in the previous analysis, F^/a is also determined in bins oiQ^ in the region 0.1 < 
X < 0.6. The new results, together with the result of the previous analysis 0, are shown 
as a function of in Fig. ^ and listed in Table ^ For the results based on the data 
taken at ^/s^ = 161 — 172 GeV the systematic error is evaluated as described above. 
The results of the two analyses at very different centre-of-mass energies, using different 
detector components to detect the scattered electron for similar values, are found to 
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be consistent. An increase of F2 /a with is observed in the data, in agreement with 
the QCD predictions and with results obtained by other experiments [0]. The data 



are compared to the LO predictions of the GRV and the SaSlD parametrisations, both 
including the contribution to from massive charm quarks, and to a higher order 
(HO) calculation |23| based on the HO GRV parametrisation for three light quarks, 
complemented by the contribution of charm quarks to F2 based on the HO calculation 
using massive charm quarks of [0], as described in more detail in The differences 



between the three predictions are small compared to the experimental errors, and all 
predictions nicely agree with the data. 

In addition the data are compared to the asymptotic prediction as detailed above. 
This approximation lies higher than the data at low and approaches the data at 
the highest reached. 

The evolution of F2 with is measured by fitting a linear function of the form 
a + h ln(Q^/GeV^) to the four data points obtained at ^/s^ = 161 — 172 GeV. Here a 
and b are parameters which do not depend on x. The result is 

F^{Q^)/a = (0.24 ± 0.11+°;?^) + (0.06 ± 0.041°;°^) In (QVGeV^) . 

The central values and statistical errors of the two parameters are obtained by a fit to 
the central values of the individual points listed in Table |^ using their statistical er- 
rors. The quality of the fit is satisfactory with x^/dof = 0.33 for the central value. The 
systematic errors are obtained from the maximum deviation from the central values to 
either side observed when repeating the fit for all variations of cuts, structure functions 
assumed in the Monte Carlo and Monte Carlo models as discussed above. This proce- 
dure automatically takes into account the correlation of the systematic uncertainties 
between the individual points. 

This measurement is consistent with the result of the previous analysis which 
used the data sample at ^/s^ = 91 GeV in the range 7.5—135 GeV^, and which 
was based on a slightly different fitting procedure. A combined fit to the OPAL data 
at ^ = 91 GeV and 161 - 172 GeV in the range of 7.5-135 GeV^ yields: 

F^{Q^)/a = (0.16 ± 0.05l°;i^) + (0.10 ± 0.021°;°^) In (gVGeV^) 

with x^/dof = 0.77 for the central value. The slope a~^dF2 /dlnQ"^ is significantly 
different from zero and the error is somewhat reduced compared to due to the 
inclusion of the data at higher centre-of-mass energies. 

The QCD prediction for the scaling violation behaviour of is different from the 
proton case. With increasing Q^, the proton structure function, F^, increases at small 
X and decreases at large x, (x > 0.1). In contrast, due to the pointlike coupling of 
the photon to quarks, the photon structure function F^ is expected to increase with 

for all values of x, and the size of the scaling violation is expected to depend on 



X. To examine whether the data exhibit the predicted variation in a ^dF2 /dlnQ"^ 



the Q range 1.86—135 GeV is analysed using common x ranges, as listed in Table 
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The published resuhs are rebinned in x but not reanalysed, i.e. the systematic error is 
evaluated as described in 

Figure ^ shows the measurement in comparison to the HO calculation. The points 
of inflection of F2 for below 15 GeV^ are caused by the charm threshold. The 
predictions are in agreement with the observed evolution of F2 with in all ranges 
of X. In order to observe experimentally the variation of a~^dF^/d InQ^ with x more 
data and a reduction of the systematic error are needed. 



7 Summary and conclusions 

The complete data samples taken by the OPAL experiment at LEP at ^/s^ = 91 GeV 
and ^/s^ = 161 — 172 GeV have been used to study the hadronic structure of the 
photon. 

A slightly better agreement between the predictions for the hadronic energy flow 
of the various models and the data in the region 0.1 < a;vis < 0.6 is found for the data 
taken at ^/s^ = 161 — 172 GeV than for the data collected at ^/s^ = 91 GeV. At 
Xyis < 0.1 significant differences persist, as the data prefer a more pointlike hadronic 
energy flow than assumed in either HERWIG or PYTHIA, and at Xvis > 0.6 the 
statistical precision is insufficient to draw firm conclusions. 

Using the OPAL data at ^/s^ = 161 — 172 GeV the photon structure function 
F2{x, Q^) has been unfolded as a function of x in four intervals, with mean momen- 
tum transfers (Q^) = 9, 14.5, 30 and 59 GeV^. Combining the OPAL data at ^/s^ = 
91 GeV and 161-172 GeV, the evolution of with in the range 0.10 < a; < 0.60 has 
been measured to be F^{Q'^)/a = (0.16 ± 0.05l|];}^) + (0.10 ± 0.02^°:°^) In (QVGeV^), 
where the first error is statistical and the second systematic. From the OPAL data 
alone the slope ^"^dFj^/dlnQ^ is found to be significantly different from zero and 
consistent with the logarithmic evolution of with expected from QCD. The de- 
pendence of the scaling violation of F2 on x was studied in three ranges in x, 0.02—0.10, 
0.10 — 0.25 and 0.25 — 0.60. The QCD prediction was found to be in agreement with 
the data, but the accuracy of the data is insufficient to show a significantly different 
slope of F2 for the three ranges in x studied. The precision of the measurement of the 
logarithmic slope of is primarily limited by systematic errors due to the modelling 
of the hadronic final state, not by the experimental statistics. There would be scope for 
considerable improvements if the Monte Carlo models could be improved to represent 
all aspects of the data. 

The data, over the x and range studied, are equally well described by sev- 
eral of the available parton density parametrisations, including the GRV and SaSlD 
parametrisations used in this analysis. They are also satisfactorily described by an F2 
based on the asymptotic solution in leading order for three light fiavours augmented 
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by the Bethe-Heitler contribution of charm quarks and a leading order hadronic part 
based on VDM. 
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Q2 range [GeV^] 


(g2) [GeV2] 


Ax 




6-11 


9 


0.020 - 


0.100 


0.33 ±0.03i|j;g^ 






0.100- 
0.250 - 


0.250 
0.600 


0.29 ±0.041°:°^ 

0.39 ±o.o8lE5:?[5 


11-20 


14.5 


0.020 - 


0.100 


0.37 ±0.03l[]:^? 






0.100 - 


0.250 


0.42 ± 0.05+o;54 






0.250 - 


0.600 


0.39 ± 0.061°:}? 


20-40 


30 


0.050 - 


0.100 


0.32 ± 0.04lo;o2 






0.100- 
0.350 - 


0.350 
0.600 


0.52 ± 0.05 1^;?^ 
0.41 ± 0.09l:o:o5 






0.600 - 


0.800 


0.46 ±0.15lgil 


40-100 


59 


0.050 - 


0.100 


0.37 ± 0.06Io;o7 






0.100- 


0.350 


0.44 ±0.071^:^? 






0.350 - 


0.600 


0.48 ±0.09tg:}^ 






0.600 - 


0.800 


0.51 ±0.14^°:^^ 


6-20 


11 


0.020 - 


0.100 


0.34 ±0.02l°:J^ 






0.100- 


0.250 


0.36 ±0.03iHo 






0.250 - 


0.600 


0.39 ±0.04^°:^^ 


20-100 


41 


0.050 - 


0.100 


0.33 ±0.03i^:J^ 






0.100- 


0.350 


0.49 ±0.04^°:°^ 






0.350 - 
0.600 - 


0.600 
0.800 


0.43 ±0.07i°:J^ 

0.54 ±o.iii^:n 



Table 1: Results for as a function of x for four active flavours in bins of for the 
data taken at ^/s^ — 161 — 172 GeV. The first errors are statistical and the second 
systematic. The SW sample with {Q"^) — 11 GeV^ consists of the data from (Q^) = 9 
and 14.5 GeV^ and the FD sample with (Q^) = 41 GeV^ consists of the data from 
(g2) = 30 and 59 GeV^. 
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Ref. 


7.5 


0.36 ±0.02 




14.7 

135 


0.41 ±0.01 
0.71 ±0.07 


t^ol i 
^-^ot i 


9 
14.5 


0.36 ± 0.05 
0.41 ±0.04 


+ 0.08 
-0.06 
+ 0.04 
-0.11 


30 

59 


0.48 ±0.05 
0.46 ±0.06 


+ 0.06 
-0.07 
+ 0.07 
-0.04 


11 
41 


0.38 ±0.03 
0.47 ±0.04 


+ 0.06 
-0.03 
+ 0.06 
-0.03 



Table 2: Results for for four active flavours in bins of for 0.10 < x < 0.60. The 
first errors are statistical and the second systematic. The bins in for the analysis 
of the data taken at y^s^ = 161 — 172 GeV are defined in Table |l[ The SW sample 
with (Q2) = 11 GeV^ consists of the data from (Q^) = g and 14.5 GeV^ and the FD 
sample with (Q^) = 41 GeV^ consists of the data from (Q^) = 30 and 59 GeV^. The 
combined results at (Q^) = 11 and 41 GeV^ are shown for completeness and are not 
used in the fit of the slope of F2, Fig. |a. 



Ax 


0.02-0.10 


0.10-0.25 


0.25-0.60 














F^Ja 




Ref. 


1.86 


0.22 ±0.01 


+ 0.06 
-0.04 










i 


3.76 


0.30 ±0.01 


+ 0.06 
-0.03 


0.32 ±0.02 


+ 0.10 

-0.03 






i 


7.5 


0.27 ±0.02 


+ 0.02 
-0.06 


0.32 ±0.02 


+ 0.08 
-0.14 


0.38 ±0.04 


+ 0.05 
-0.18 


i 


14.7 


0.37 ±0.01 


+ 0.05 
-0.15 


0.40 ±0.02 


+ 0.06 
-0.04 


0.42 ± 0.02 


+ 0.09 
-0.09 


i 


135 






0.62 ±0.09 


+ 0.12 

-0.05 


0.73 ±0.07 


+ 0.04 
-0.05 


i 


11 


0.34 ±0.02 


+ 0.11 
-0.02 


0.36 ±0.03 


+ 0.02 
-0.10 


0.39 ±0.04 


+ 0.12 
-0.04 




41 


0.32 ±0.03 


+ 0.19 
-0.03 


0.46 ±0.04 


+ 0.02 
-0.13 


0.46 ±0.05 


+ 0.10 
-0.06 





Table 3: Results for for four active flavours in bins of and x. The first errors are 
statistical and the second systematic. The bins in for the analysis of the data taken 
at ^/s^ =161 — 172 GeV are defined in Table |l|. The results from are rebinned 
to match the same ranges in x. For statistical reasons only the combined samples with 
(Q^) = 11 and 41 GeV^ are used for the subdivision into several bins in x, Fig. 
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Figure 1: The hadronic energy flow per event as a function of pseudorapidity t] for 
the data taken at ^/s^ = 161 — 172 GeV and various signal Monte Carlo samples. The 
energy flow is shown for two ranges of Xvis for the SW and FD samples. The errors 
shown are statistical only. The tagged electron is always at negative pseudorapidity 
and is not shown. The different regions in rapidity indicated denote the forward region 
(FR) which contains the SW and FD detectors and the central detector (CD). 



19 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



r 



< = 11 GeV^ 

-i- OPAL data 



— Y*Y"^h*drons unf. 
* + . 

^ ^ 

IZ3 y y^e e 




(a) 

■ ' ■ ' I I i_ 



0.2 0.4 0.6 0.8 



j2 100 F 
I 90 j- 
80 I 

Y*Y-^hadrons J 70 

60 

50 ^ 
40 
30 

20 ^ 
10 % 




T — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



< Q"> = 41 GeV 

OPAL (b) 



1 — I — rq 
2 




0.2 0.4 0.6 0.8 



T — I — I — I — I — I — I — I — I — I — I — I — I — r 







< Q^> = 11 GeV^ 

OPAL (c) 




g 100 



20 30 

[GeV^] 



80 
60 
40 
20 




T — I — I — I — I — I — I — I — I — I — I — I — I — r 







< Q^> = 41 GeV^ 

OPAL (d) 




100 



150 

[GeV^] 



Figure 2: The Xyis, (a) and (b), and Q^, (c) and (d), distribution of the SW and 

FD samples for the data taken at y^s^ = 161 — 172 GeV. The dashed histogram, 
7*7 hadrons, shows the events from the HERWIG Monte Carlo, using the GRV 
parametrisation and the standard cuts (the reference sample), with the background 
added, before the unfolding; the solid histogram, 7*7 — > hadrons unf., shows the same 
quantity, but reweighted based on the result of the unfolding. The background events 
from the reactions 7*7 t^t~ and 7*7 e^e~ are also shown separately at the 
bottom of the figure. All errors shown are statistical only. 
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Figure 3: The measurement of for the data taken at ^/s^ — 161 — 172 GeV. The 
structure function is measured for four active flavours in four bins in with mean 
values of (a) (Q^) = 9 QeV^, (b) (Q^) = 30 GeV^, (c) (Q^) = 14.5 GcV^, and (d) {Q^) 
= 59 GeV^. In (e) the measurement for the combined data sets of (a) and (c), and in (f) 
the measurement for the combined data sets of (b) and (d) is shown. The points show 
the measured F^. The bin sizes are indicated by the vertical lines at the top of the 
figure. The solid line represents the derived from the GRV parametrisation and the 
dot-dashed line denotes the F2 derived from the SaSlD parametrisation, both using the 
Bethe-Heitler contribution to F2 for massive charm quarks. The charm quark mass rUc 
is taken to be 1.3 GeV and 1.5 GeV in the case of SaSlD and GRV, respectively. The 
dotted fine (HAD) represents the hadronic component and the dashed line (ASYM) 
the full augmented asymptotic F^, as described in the text. 
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Figure 4: The measurement of for four active flavours as a function of Q^, (a) for 
the range 0.1 < a; < 0.6, and (b) subdivided into 0.02 < x < 0.10, 0.10 < a; < 0.25 
and 0.25 < x < 0.60. In addition shown in (a) are the of the GRV (LO) and the 
SaSlD (LO) parametrisation, the of the asymptotic prediction (ASYM) and the 
result of a higher order calculation (HO), where the last two predictions are only shown 
for > 4 GeV^. In (b) the data are only compared to the HO prediction. In both 
figures the errors are statistical and systematic. In some of the cases the statistical 
errors are not visible because they are smaller than the size of the symbols. 
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